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ABSTRACT 

We present a comparison of optical and X-ray properties of galaxy clusters in the northern 
sky, using literature data from BAX and optically selected clusters in DPOSS. We determine 
the recovery rate of X-ray detected clusters in the optical as a function of richness, redshift 
and X-ray luminosity, showing that the missed clusters are typically low contrast systems when 
observed optically (either poor or at high redshifts). We employ four different statistical tests 
to test for the presence of substructure using optical two-dimensional data. We find that ap- 
proximately 35% of the clusters show strong signs of substructure in the optical. However, the 
results are test-dependent, with variations also due to the magnitude range and radius utilized. 
We have also performed a comparison of X-ray luminosity and temperature with optical galaxy 
counts (richness). We find that the slope and scatter of the relations between richness and the 
X-ray properties are heavily dependent on the density contrast of the clusters. The selection of 
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substnicturc-frcc systems does not improve the correlation between X-ray luminosity and rich- 
ness, but this comparison also shows much larger scatter than one obtained using the X-ray 
temperature. In the latter case, the sample is significantly reduced because temperature mea- 
surements are available only for the most massive (and thus high contrast) systems. However, 
the comparison between temperature and richness is very sensitive to the exclusion of clusters 
showing signs of substructure. The correlation of X-ray luminosity and richness is based on the 
largest sample to date 750 clusters), while tests involving temperature use a similar number 
of objects as previous works ( ^ 100). The results presented here are in good agreement with 
existing literature. 

Subject headings: galaxies: clusters: general - X-rays: gala^xies - catalogs 



1. Introduction 

In most scenarios of structure formation, galaxy 
clusters are the largest and latest objects to form 
under the influence of their own gravity. Reliable 
measurements of cluster masses allow us to study 
the evolution of the large scale structure through 
the cluster mass function and its time evolution. 
Several studies have examined the evolution ex- 
perienced by clusters, using the results to provide 
estimates of important cosmological parameters, 
such as ftm and cts (Eke et al. 1998; Mathicscn and 
Evrard 1998; Kitayama and Suto 1997; Bahcall et 
al. 1997, 2003; Carlberg et al. 1997; Donahue et al. 
1998: Rcichart et al. 1999; Blanchard and Bartlctt 
1998; Postman, Lauer, Oegerle, & Donahue 2002). 
The main difficulty faced by these studies lies in 
determining the cluster masses. A variety of tech- 
niques can be used for this measurement, including 
dynamical methods (velocity dispersions), mea- 
suring the temperature of the intra-cluster gas, 
observations of the Sunyaev— Zeldovich effect, and 
weak gravitational lensing. All of these methods 
arc observationally expensive, which render large 
scale structure studies based on extensive samples 
impractical. 

An alternative method to obtain mass estimates 
for a large sample of clusters is to connect the clus- 
ter mass to an easily observable parameter, such 
as X-ray luminosity (Lx) or the number of con- 
stituent galaxies (richness). The tightest mass- 
observable relations involve X-ray temperatures 
(Tx) and spectroscopically measured velocity dis- 
persions {ay). Unfortunately, large samples of 
clusters with either or both of these parameters 
measured do not exist currently. The X-ray lu- 
minosity (Lx) and optical richness (or luminos- 
ity, Lopt) also correlate well with mass, although 
with larger scatter than relations with T^ and cr„. 



The precise; ni(;a.surement and calibration of these 
relations are essential for reliable mass function 
estimates. For instance, underestimation of the 
scatter in the mass— observable relation can result 
in overestimation of erg (Voit 2005). Furthermore, 
the correlation of mass to properties determined 
at different wavelengths can lead to contradictory 
results, which may be associated with the selec- 
tion function of the cluster catalogs constructed 
in different regimes, or with physical processes and 
evolution within the cluster populations. 

Thus, understanding the systematics present in 
optical and X-ray surveys, and the way they com- 
plement each other, are necessary precursors to 
studies of cluster evolution and are even more im- 
portant for cosmological tests that require statisti- 
cally robust cluster samples. Some of these issues 
have been addressed in recent literature, either by 
comparing X-ray and optical cluster catalogs or 
conducting joint X-ray/optical surveys of galaxy 
clusters (Donahue et al. 2001, 2002; Gilbank 2001; 
Gilbank et al. 2004; Yee and EUingson 2003; Basi- 
lakos et al. 2004; Popesso et al. 2004, 2005; Smith 
et al. 2005). A very important by-product of such 
comparisons is the selection and follow-up study 
of unusual clusters, such as those that are opti- 
cally rich but X-ray underluminous (Gilbank 2001; 
Gilbank et al. 2004; Lubin et al. 2004). 

As mentioned above, an important issue con- 
cerning these scaling relations is the precision to 
which they can be established and which observa- 
tional biases contribute to increasing their scatter. 
Information on the distribution of galaxies, hot 
gas and dark matter within clusters can play an 
important role in understanding this scatter. The 
diversity in the dynamical states of galaxy clusters 
may be associated with the uncertainties in the de- 
termination of the scaling relations. The presence 
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of substructure is a clear sign of incomplete relax- 
ation in a cluster, and can be estimated from the 
X-ray emission from the intra-cluster medium as 
well as the distribution (one, two or three dimen- 
sional) of cluster galaxies (Mathiesen, Evrard and 
Mohr 1999; Kolokotronis et al. 2001; Smith et al. 
2005). The dynamical state of galaxy clusters is 
also strongly related to the underlying cosmology, 
so that the evolution of cluster substructure with 
look-back time is a powerful tool for constraining 
cosmological parameters (Mohr et al. 1995). 

A few recent works have investigated the pos- 
sibility that cluster substructure adds scatter to 
the observed scaling relations, but the results 
have been ambiguous. While Smith et al. (2005) 
find that the observational scaling relations of 
clusters behave differently for relaxed and unre- 
laxed clusters, O'Hara et al. (2006) find the oppo- 
site (based both on simulations and observational 
data). Some of these differences may be due to the 
relatively small sample size or to the inconsistent 
techniques for measuring substructure. 

The main goal of this paper is the establishment 
of scaling relations for galaxy clusters through the 
comparison of optical and X-ray properties {l^gais, 
Lx, Tx). We also want to examine which factors 
contribute to the observed scatter in these rela- 
tions. We first investigate the recovery rate of X- 
ray emitting galaxy clusters in an optical survey, 
searching specifically for X-ray detected clusters 
which could be missed in the optical. We also 
conduct a substructure study based solely on the 
photometric data. Finally, wc compare the optical 
and X-ray properties of galaxy clusters, investigat- 
ing the dependence of the results with the contrast 
cut of the sample used, the centroid adopted for 
the richness measures and the amount of substruc- 
ture. We limit this study to photometric and X- 
ray parameters and investigate the systematics in 
these comparisons, postponing to a future work 
the extension of the cluster scaling relations for 
other parameters, such as i?200) cr^ and especially 
mass. 

We use a list of X-ray observed galaxy clus- 
ters from BAX (http://bax.ast.obs-mip.fr/) as our 
basis. The catalog initially comprises 914 galaxy 
clusters in the northern hemisphere, spanning the 
redshift range 0.05 < z < 0.40. We trim this hst to 
the same area covered by the Northern Sky Opti- 
cal Cluster Survey (NoSOCS) in order to compare 



the two catalogs. The final X-ray list contains 792 
clusters, of which 638 (81%) are detected in the 
optical. 

The remainder of this paper is organized as fol- 
lows. In §2 we describe the X-ray data, while in 
§3 we provide details about the NoSOCS cata- 
log. The BAX and NoSOCS comparison is dis- 
cussed in §4, while the substructure analysis is 
presented in §5. In §6 we investigate the correla- 
tion between optical and X-ray properties, namely 
Ngais, Ljjc and Tx- The X-ray temperature- 
luminosity relation is also determined. In particu- 
lar we investigate the impact of substructure and 
the cluster contrast on the scatter of these rela- 
tions. We summarize our results in §7. We use 
Ho = 100 h km s^^ Mpc~-^, except where the 
/i— dependence is explicitly given, and qq = 0.5 {h 
is fixed to 0.5 only for converting magnitudes, but 
not for computing angular sizes). In order to be 
in agreement with BAX (see http://bax.ast.obs- 
mip.fr/litnil/lielp/baxJielp.litml for notes on the 
BAX data) we decided not to adopt the currently 
favored cosmology with h = 0.70, flm = 0.30 and 
Qa =0.70. This choice has no significant effect on 
the results for the scaling relations (namely slope, 
intercept and scatter). The only parameters that 
might be affected by changing cosmology are those 
associated with evolution, which is not expected to 
be significant in the current work due to the small 
redshift range probed. 

2. X-ray Data 

We search for optical counterparts of X-ray 
observed clusters selected from BAX, Base de 
Donnees Amas de Galaxies X (http://bax.ast.obs- 
mip.fr/). BAX is an on line research database con- 
taining information on all galaxy clusters with X- 
ray observations to date. BAX contains all galaxy 
clusters confirmed as X-ray sources, but is not re- 
stricted to objects first detected in X-rays, and 
includes optically-selected clusters (such as Abell 
clusters) with X-ray follow-up. The database 
allows us to retrieve each cluster's coordinates 
and redshift; various X-ray observational measure- 
ments, namely flux (Fx), luminosity (Lx) and 
temperature (Tx)', and a set of bibliographical ref- 
erences. The cluster coordinates and redshift are 
automatically generated from the NASA/IPAC 
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Extragalactic Database (NED^). Fx and Lx are 
in the ROSAT [0.1-2.4 keV] band. Fx and Lx 
units are 10~^^ ergs cm~^ and 10^^ ergs s~^, 
respectively. Only BAX canonical values of Fx, 
Lx and Tx are retrieved. These are values chosen 
among the most precise and/or latest published 
measurements, meaning that if there are miilti- 
ple references for a given cluster, the most precise 
and/or recent is kept. The BAX output is given 
for Ho = 50 km s"^ Mpc'^ and f7„ = 1.0. At 
the time of our query, BAX contained information 
on 1851 groups and clusters of galaxies, with 1656 
clusters having luminosity measurements and 463 
clusters with temperature measurements. 

We performed a multi-criterion search with the 
following constraints: —5° < ^ < 90° and 0.05 < 
z < 0.40. These limits ensure overlap with our 
optical galaxy cluster catalog (NoSOCS), which is 
based on DPOSS {S > -2.5° and z < 0.3). The 
query was performed on 2005 April 12, resulting 
in a list of 914 galaxy clusters. We then com- 
pared the X-ray centers of these clusters to all the 
plate limits used for NoSOCS, keeping only those 
that fall within these limits. We found 800 X-ray 
clusters within the NoSOCS limits, from which we 
eliminated eight that overlapped with bad areas 
from our survey (excised regions due to saturated 
objects such as bright stars). Of these 792 re- 
maining systems, 638 have an optical counterpart 
in our catalog. From these 638 common clusters, 
620 have measured X-ray luminosities (with two 
others discarded due to problems measuring the 
optical richness), while 101 have X-ray tempera- 
tures determined. More details on the recovery 
rate of X-ray galaxy clusters by DPOSS are pro- 
vided in section §4. The comparison of optical and 
X-ray properties, shown in section §6, is based on 
this list of 618 common clusters with available Lx 
and Nga/s, or the subset of 101 clusters for which 
temperatures are available. 

3. Optical Data 

The X-ray galaxy cluster catalog derived from 
BAX is compared to the Northern Sky Optical 
Cluster Survey (NoSOCS, Gal et al. 2003, 2006; 
papers II and III, respectively), which is a cata- 
log of galaxy clusters constructed from the digi- 
tized version of the Second Palomar Observatory 
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Sky Survey (POSS-II; DPOSS). NoSOCS is de- 
rived from high-latitude fields |6| > 30°, cover- 
ing ~ 11000 deg^ and containing ~ 16000 cluster 
candidates. Richness and redshift estimates are 
provided for nearly all clusters; the richness mea- 
sure (denoted l^gais) is described in detail in Gal 
et al. (2003) and in the next section. The me- 
dian redshift of the sample is Zmed = 0.16, with 
a median richness of l^gais,med = 31. Prom the 
selection functions presented in papers II and III 
we expect to find rich clusters (Ngo/g ~ 100) out 
to z 0.3. This is the largest galaxy cluster cat- 
alog available to the present date. The NoSOCS 
catalog and other DPOSS products can be found 
at http://dposs.ncsa.uiuc.edu/. 

4. Richness and the BAX - NoSOCS 
Comparison 

Due to differences between the photometric red- 
shift estimates from NoSOCS and the spectro- 
scopic redshifts provided by BAX, we adopt the 
BAX redshift for all analysis involving clusters 
common to the two datasets. We note that the 
scatter in the comparison of NoSOCS and BAX 
redshifts is comparable to the redshift accuracy es- 
timated by Gal et al. (2003). We utilize the BAX 
redshifts measures as they are generally more ac- 
curate (as they are typically spectroscopic), and 
the richness measure and its accuracy are sensi- 
tive to the cluster centroid, radius and redshift 
adopted. This choice is corroborated by the fact 
that the comparisons of optical and X-ray prop- 
erties show reduced scatter when using the BAX 
redshifts (described in Sections 4.1 and 6). The 
improved correlations are expected, since we now 
compute optical richnesses with the same redshifts 
used for calculating the X-ray properties (Lx and 
Tx taken directly from BAX). In §4.1 we de- 
tail the richness estimator, investigating its depen- 
dence on the cluster centroid, redshift and aper- 
ture used. Then, in §4.2 we show the results - 
in terms of recovery rate - of the comparison be- 
tween the BAX clusters and the NoSOCS. We be- 
gin by searching for possible optical counterparts 
in NoSOCS for each of the 792 BAX clusters, using 
a matching radius of 1.50 h~^ Mpc and keeping 
only the closest match to each BAX entry. The 
total number of recovered systems is 638. 
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Fig. 1. — Richness residual for different cluster centroids and adopted redshifts. The labels in the lower 
right of each figure have the following meaning: "C" stands for coordinates and "Z" for redshifts; "O" is the 
original cluster coordinate; "B" is the BAX ccntroid or rcdshift; "L" is the luminosity weighted coordinate; 
"P" is the centroid given by the maximum density peak; and "D" means the photometric redshift from 
DPOSS. Richness is compared for the following cases: (a) cluster original center vs. BAX centroid, both 
using the BAX rcdshift: (b) using the DPOSS redshift vs the BAX rcdshift (both at the original coordinates); 
(c) original center vs the luminosity weighted coordinate (both using the BAX redshift); (d) original center 
vs the density peak centroid (both using the BAX redshift); (e) same as panel (c), but for the 10190 clusters 
used to estimate substructure; (/) same as panel (d), but for all 10190 clusters in the NoSOCS sample. Panels 
(e) and (f) use the DPOSS photometric redshift by necessity. 



4.1. Systematic Effects on Richness Esti- 
mation 

Measurement of the optical richness of a cluster 
depends on the cluster centroid, radius adopted, 
and the redshift. Our richness measure Ngau is 
defined as the number of galaxies at m* — 1 < 
nir < m* + 2 within a given aperture, where m* is 



the characteristic apparent magnitude of the clus- 
ter luminosity function. For this work we assume 
a universal Schecter luminosity function, with pa- 
rameters given by Paolillo et al. (2001). The 
characteristic magnitude is M* = —22.17, while 
a = — 1.1. As explained below, we find the opti- 
mal aperture to be 0.50 h~^ Mpc. We have consid- 
ered no evolution in the characteristic magnitude 
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(m*). k— corrections are applied to the magnitude 
of each galaxy, as explained in items (1) and (3) 
below. These corrections are obtained trhough the 
convolution of spectral energy distributions from 
Coleman et al. (1980) with the DPOSS r filter. 
The steps taken to estimate richness are: 

1. We use the redshift of each cluster to con- 
vert the characteristic magnitude to an ap- 
parent magnitude (m*) and to calculate the 
apparent radius (in seconds of arc) for a fixed 
aperture of 0.50 h~^ Mpc. If a cluster lies 
closer than 0.50 h~^ Mpc to a plate border 
we do not calculate its richness. The redshift 
is also used to compute the k— correction val- 
ues typical of elliptical and late-type galax- 
ies (Sbc) at the cluster redshift. These are 
named "ke" and "ks", respectively. Be- 
cause we only want to count galaxies with 
m* — 1 < rrir < m* + 2, we select all galax- 
ies within 0.50 h~^ Mpc of the cluster cen- 
ter at m* — 1 + ks < rrir < m* + 2 + ke. 
The k— corrections are applied to individual 
galaxies at a later stage, so these limits guar- 
antee that we select all galaxies that can fall 
within m* — 1 < m.r < m* + 2. The num- 
ber of galaxies selected in the cluster region 
is Ncfu. If the low (m* — 1 + ks) or high 
(m* + 2 + ke) magnitude limits fall outside 
the survey observational limits (15.0 < rur < 
20.0), we apply a correction to the richness 
measure (explained below). 

2. We estimate the background contribution 
locally. We randomly select twenty back- 
ground boxes (1200" X 1200") within an an- 
nulus ranging from 2.25 h~^ Mpc to 2.25 h~^ 
Mpc -|- 1.3° of the cluster center. Galax- 
ies are selected within the same magnitude 
range as used for computing Nc;„. We re- 
select boxes that overlap with an excised 
area of the survey (due to bright stars, for in- 
stance). The background counts of each box 
are scaled to the same area as the cluster, 
and the median counts (from all 20 boxes) 
gives the background estimate {N^kg). We 
adopt the interquartile range (IQR, which is 
the range between the first and third quar- 
tiles) as a measure of the error in l>ibkg , which 
we term Qa^kg (foi' normally distributed data 
IQR = 1.35 XCT, where a is the standard de- 



viation). The background corrected cluster 
counts (Nciu - Nbfeg) is called Ncorr- 

3. Next, a bootstrap procedure is used to sta- 
tistically apply k— corrections to the galaxy 
populations in each cluster. In each of 100 
iterations, we randomly select Ncorr galax- 
ies from those falling in the cluster region 
(Nc/u). We then apply a k— correction to 
the magnitude of each galaxy. An elliptical 
k— correction is applied to 80% of the Ncorr 
galaxies, while an Sbc k— correction is ap- 
plied to the remaining 20%. Finally, we use 
these k— corrected magnitudes to count the 
number of galaxies at m* — l < rrir ^ m* + 2. 
The final richness estimate Ngais is given by 
the median counts from the 100 iterations. 
The richness error from the bootstrap pro- 
cedure alone is given by Quboot — IQR- The 
total richness error is the combination of this 
error and the background contribution, so 

that = ./oiZ+QlZ- 

4. If the cluster is too nearby or too distant, 
either the bright (m* — 1 -|- ks) or faint 
(m* -1-2-1- ke) magnitude limit, respectively, 
will exceed one of the survey limits (15.0 < 
rUr < 20.0). In this case we need to apply 
one of the following correction factors to the 
richness estimate: 

, , O m dm 

Jm* — 1 ^ ' ^..N 

~ rm'+2 -r, \ / 

Jj^g 9(ni)am 

Xn.-i '^{m)dm 
T2 = 720 ^, (2) 

We call 71 and 72 the low and high magni- 
tude limit correction factors. Whenever nec- 
essary, one of the above factors is multiplied 
by Nga(s. Because the optical data span 
five magnitudes it is impossible to simultane- 
ously exceed both the bright and faint end. 

An estimate of the cluster contrast is also ob- 
tained along with the richness. The contrast is 
defined as the ratio between the number of galax- 
ies in the cluster region having m* — 1 + ks < rrir < 
m*+2+ke {^du), and the background error Qa^kg 
(estimated within the same magnitude range). 
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Wc utilize this richness measure to examine its 
dependence on the cluster centroid, radius and 
redshift. In Figure 1 we show the residual be- 
tween different richness estimates obtained when 
adopting different centroids and redshifts. In the 
lower right corner of each panel, we encode the 
coordinates "C" and redshifts "Z" used as follows: 

• CO: the original cluster position from NoSOCS, 
which is the location provided by SExtractor 
in an adaptively smoothed density map; 

• CL: a luminosity weighted coordinate, where 
each galaxy position within the cluster aper- 
ture is weighted by the galaxy's luminosity 
and a new center is computed; 

• CB: the cluster center cataloged in BAX; 

• CP: the coordinate of the max;imum density 

peak; 

• ZD: the photometric redshift from the 
DPOSS data; 

• ZB: the redshift (usually but not always 
spectroscopic) cataloged in BAX. 

In the two lower panels we show the results ob- 
tained with all 10190 clusters used to estimate 
substructure (section 5), while in the upper four 
panels the analysis is restricted to the 636 clusters 
common to BAX and DPOSS. 

Prom inspection of this figure it is clear that the 
richnesses estimated with the X-ray centroids are 
underestimated compared to those obtained with 
the optical centers (panel a). This is likely due 
to the optical center being a better choice than 
the X-ray center to trace the distribution of galax- 
ies within the cluster. If we restrict the study to 
clusters with small optical-X-ray centroid offsets 
( < 0.50 h"^ Mpc) the results are very similar, 
although a small trend for higher richness val- 
ues computed with the optical centers is expected. 
When adopting the same center (optical) and us- 
ing different redshifts (from BAX and DPOSS) 
there is no obvious trend (the median residual 
is zero), but the scatter is the largest of all the 
comparisons (panel b). This implies that there is 
no systematic difference between the photometric 
redshifts and those from BAX. The scatter arises 
from a shift in sampling the luminosity function 
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Fig. 2. — Top: dependence of the linear correlation 
coefficient (p) for the the hx—^gais (solid line) 
and Tjsf— Nga/s relations (dotted line) as a function 
of the radius; Bottom: variation of the orthogonal 
scatter of the \jx—^gais (solid line) and Tx—^gais 
relations (dotted line) as a function of the radius 
used for calculating richness. 

of each cluster and a change in radius when esti- 
mating richness. Adopting either the original or 
the luminosity weighted coordinates does not ap- 
preciably affect the resultant richnesses (panels c 
and e), as the two centers are typically similar. 
On the other hand, the richness estimate obtained 
with the density peak center is underestimated in 
comparison to the other optical centroids (panels 
d and f). This may be due to the presence of sub- 
structure (or projection effects), as a cluster show- 
ing more than one galaxy concentration within the 
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Table 1 

Scatter and number of objects in the optical vs X-ray relations for different richness 

estimates 



Set 


Lj!f-Nga(s 










All clusters 


/?sigm/>0.06 


All clusters 


l3signif>0.06 


CO-ZB 


0.172; 430 


0.171; 271 


0.112; 53 


0.041; 23 


CB-ZB 


0.184; 418 


0.174; 233 


0.119; 53 


0.061; 26 


CO-ZD 


0.179; 424 


0.175; 266 


0.111; 48 


0.111; 29 


CL-ZB 


0.174; 432 


0.174; 274 


0.120; .53 


0.062; 28 


CI'-ZB 


0.171; 120 


0,17."): 2()0 


0.120; 5:-! 


0.0()5: 20 



measurement aperture (0.50 Mpc) will have 
multiple density peaks. Choosing the highest den- 
sity peak for richness calculation i.e., counting 
galaxies within 0.50 h~^ Mpc from this new cen- 
ter - may not be appropriate for tracing the galaxy 
distribution. If the cluster is defined by a single 
peak, then the richness estimate should not be af- 
fected. 

Beyond considering the effects mentioned above, 
we investigate which metric aperture is optimal for 
computing the richness. In Gal et al. (2003) we 



{N,.,.-,„„ / N,.J vs. N,.,. 
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Fig. 3. — Relative richness error as a function 
of richness for the full NoSOCS sample (~ 16500 
clusters) . The solid line is the variation of relative 
richness error using a 0.50 Mpc aperture for 
computing richness. The variation when adopting 
a radius of 1.00 Mpc is shown by the dotted 
line. 



adopted a 1.00 h~^ Mpc radius, which may not 
be the best choice for the current purposes. Like 
Popcsso ct al. (2004) wc used the scaling relations 
between optical and X-ray properties to guide 
the choice of the optimal radius. We investigated 
how the scatter of these relations is affected by 
the different radii used for computing richness, 
examining the Lx—^gata and Tx—^gata relations 
based on all clusters common to BAX and DPOSS 
used in the substructure analysis (see section 5). 
Figure 2 shows the variation of the linear corre- 
lation coefficient (p) and the orthogonal scatter 
(cr) for these two relations with the aperture used 
to estimate richness. We see a clear trend for 



Offset (arcsec) 
400 600 




E 100 - 



Offset (h-" Mpc) 



Fig. 4. — Offset distribution between BAX and 
NoSOCS clusters in h~^ Mpc (solid line) and arc- 
seconds (dotted line). 
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minimum scatter and highest correlation for an 
aperture of ~ 0.50 h~^ Mpc. In Figure 3 we show 
the relation between the relative richness error 
gals-err/^ gals) and richness for two apertures. 
We see that clusters with ISigais < 10 (within 0.50 
h~^ Mpc) have errors comparable to their rich- 
nesses. The increase in the richness error as we 
increase the radius from 0.50 h~^ Mpc to 1.00 h~^ 
Mpc is also striking. Based on these three plots 
and on Table 1 (where we show the scatter of the 
optical vs X-ray relations for richness estimates 
based on different centroids and redshifts) we re- 
compute the richness estimate (Ng,j;,,), instead of 
using the values from Papers II and III. We adopt 
the original optical cluster position, the BAX red- 
shift and an aperture of 0.50 h"'^ Mpc for each 
cluster, and use the BAX cosmology. We also ex- 
clude clusters with 'Ngais < 10 (unless otherwise 
stated) for the comparisons shown in §6. This step 
(re-computing richness) is crucial for minimizing 
the scatter in the comparisons between optical and 
X-ray properties (§6 and Table 1). The impact of 
centroid, radius and redshift on substructure mea- 
surement is discussed in §5. 

4.2. Comparison of the two catalogs 

We compare the BAX and NoSOCS catalogs 
using BAX as the reference, searching for opti- 
cal counterparts in NoSOCS. We adopt a search 
radius of 1.50 h~^ Mpc, keeping only the near- 
est match to each of the 792 BAX clusters. The 
total number of recovered systems is 638. In Fig- 
ure 4 we show the offset distribution in h~^ Mpc 
(solid line) and arcseconds (dotted line). The typ- 
ical X-ray— optical offset is < 0.50 h"-'^ Mpc (or 
< 400"). The redshift and richness distributions 
of the common and missing clusters are shown in 
Figure 5, as the top left and right panels, respec- 
tively. The majority of BAX clusters not identi- 
fied by NoSOCS are at higher redshifts {z > 0.2) 
and have low richnesses (Ngais < 40) systems. For 
clusters not found by NoSOCS, we adopt the X- 
ray coordinates (given in BAX) to compute their 
richnesses. 

The overall recovery rate of BAX clusters by 
NoSOCS is 81%. Considering the redshift and 
richness (or Lx) ranges sampled this is an encour- 
aging result. It is also important to keep in mind 
that the BAX catalog is a compilation of all X-ray 
galaxy clusters found in the literature and not a 



complete list to a given flux limit. Instead of try- 
ing to explain each missing cluster on a case-by- 
case basis, we investigate their overall properties 
(richness and Lx) 

clS cL fimction of redshift. This 
approach also enables us to directly compare our 
results to the selection functions (SF) estimated by 
Gal et al. (2003, 2006). The recovery rate of X-ray 
emitting clusters by NoSOCS is shown in the lower 
panels of Figure 5. In the left panel we show the 
detection rate for different richness classes, while 
in the right panel the relation is shown for different 
X-ray luminosities. Comparison of the left panel 
of this figure to Figure 6 from Paper II suggests 
that the earlier simulations performed to evaluate 
the SFs tend to underestimate the successful detec- 
tion rate. However, the BAX list used here is not 
a complete catalog and includes clusters that were 
first found in the optical and only later observed 
at X-ray wavelengths. Even considering these two 
points the results shown in the lower panels of Fig- 
ure 5 are quite impressive, demonstrating a recov- 
ery rate of approximately 90% of all X-ray lumi- 
nous {Lx ^ 3.2 10'*^ ergs s~^) galaxy clusters at 
any redshift out to ^; 0.2. In terms of richness, 
we recover all X-ray clusters cataloged in BAX 
that have Ngais <^ 80 out to z ~ 0.3, and more 
than 90% of those with 'Ngau > 25 to z ~ 0.2. 

5. Optical Substructure 

We have used photometric data alone to ob- 
tain an estimate of the fraction of clusters with 
evidence for substructure, based on our sam- 
ple of 10190 clusters at 0.07 < z < 0.21. In or- 
der to evaluate substructure we apply four two- 
dimensional (spatial) statistical tests (Pinkney et 
al. 1996). It is well known that any substructure 
test statistic has little meaning if not properly 
normalized, which can be achieved by compar- 
ing the results for the input data to those for 
substructure- free samples (the null hypothesis). 
For the 2D tests employed in this work the null 
hypothesis is given by an azimuthally symmetric, 
smooth distribution of points, with surface density 
decreasing as a function of radius. We generate 
such data via azimuthal randomization, where the 
distance of each galaxy to the cluster center is 
maintained and a new azimuth is randomly as- 
signed. The main advantage of this technique is 
the exact replication of the radial profile of the 
cluster. The four tests used here each provide a 



9 





Fig. 5. — The recovery rate of BAX clusters in NoSOCS. (a) redshift distribution for common clusters 
and BAX clusters not found in DPOSS; (6) the richness distributions for common objects and BAX only 
sources; (c) recovery rate of X-ray clusters by NoSOCS, as a function of redshift, for different richness classes 
(1 < Ngais < 25; 25 < Ng^is < 50; 50 < Ng^is < 80; 80 < Ng^is < 120); (d) the same as in panel (c), but for 
different X-ray luminosity ranges (—2.0 < LogLx < —0.5; —0.5 < LogLx < 0.5; 0.5 < LogLx < 2.0). 



significance level, given as the probability that the 
observed distribution is drawn randomly from one 
free of substructure. 

The significance level is determined through 
Monte Carlo simulations. For each input data 
file wc generate N — 500 simulated data sets by 
azimuthal randomization. We then calculate the 
number of Monte Carlo simulations which show 
more substructure than the real data. Finally, this 



number is divided by the number of Monte Carlo 
simulations. For most of the analysis in this pa- 
per we set our significance threshold at 5%, mean- 
ing that only 25 of the 500 simulated datasets can 
have substructure statistics higher than the obser- 
vations to consider a substructure estimate signifi- 
cant. This choice of threshold is explained in §5.5. 

In the next four subsections we give a brief de- 
scription of each substructure estimator: the an- 
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gular separation test (AST), the Fourier elonga- 
tion test (FE), the Lee statistic (Lee 2D) and the 
symmetry test (/3). Detailed descriptions of all four 
tests are provided by Pinkney et al. (1996). 

5.1. The Angular Separation Test 

The angular separation test was developed by 
West, Ocmlcr and Dckcl (1988) and considers the 
galaxy distribution in the cluster only in terms of 
angular coordinates around the cluster center. In 
this case, the presence of substructure is detected 
as an excess number of small angular separations 
between galaxy pairs relative to the expected num- 
ber in a spherically symmetric, substructure-free 
distribution. The harmonic mean angular separa- 
tion is given by 



2/7V(7V- 1)^0^1 



i>j 



(1) 



where the sum is performed over all pairs. N is the 
total number of galaxies and % is the angular sep- 
aration between galaxies i and j . The test statistic 
is the ratio of Ohm measured for the cluster and the 
harmonic mean for a Poisson distribution with the 
same number of galaxies, AST — Ohm/Qp- This ra- 
tio is near unity for substructure-free systems, and 
less than 1.0 for clumpy distributions. 6'p is ob- 
tained from Monte Carlo simulations whosc^ sets 
are generated using the azimutlial randomization 
described above. The test also adopts a small an- 
gle filter, rejecting all Ojj that arc; less than 1% of 
the expected mean inter-galaxy separation if the N 
galaxies were imiformly distributed (0.01[27r/N]). 
W(«t, Oc!nil(?r and T)vkv\ (1988) find that this test 
is too sensitive to the shape of the cluster density 
profile to provide unambiguous evidence for the 
presence of substructure. The main disadvantage 
of this test is the loss of information when going 
from two spatial dimensions to a single angular di- 
mension. However, as pointed out by Pinkney et 
al. (1996), AST may be a more sensitive diagnostic 
for clumping in multiple systems. 

5.2. The Symmetry Test 

The /3 (or symmetry) test was also introduced 
by West, Oemler and Dekel (1988) to test for sig- 
nificant deviations from mirror symmetry about 
the cluster center. This test assumes that a sub- 
cluster represents a local asymmetry superposed 



on an otherwise symmetric distribution. For each 
galaxy i, a local density estimate is obtained by 
the mean distance to the N'^/'^ nearest neighbors 
di (West, Oemler and Dekel (1988) used the five 
nearest neighbors). The local density for a point 
o diametrically opposite to galaxy i is estimated 
in the same way (estimating the mean distance to 
the iV^/^ nearest galaxies, do). For a symmetric 
galaxy distribution the values of di and do should, 
on average, be approximately equal, but they will 
differ for clumpy distributions. The asymmetry 
for a given galaxy i is given by 



(2) 



A = logw ( -j; 



The /? statistic is then defined by the average 
value <f3i> over all galaxies. For a symmetric dis- 
tribution < (3 >« 0, while values of </3> greater 
than indicate asymmetries. </?> is a density- 
weighted average, being more sensitive to the pres- 
ence of substructure than if random points were 
selected from the cluster. The /3 test is sensitive 
to deviations from mirror symmetry, but not nec- 
essarily to ones from circular symmetry. In other 
words, the test does not mistake an elongated 
smooth cluster for one containing substructure. 

5.3. The Fourier Elongation Test 

One of the characteristics of galaxy clusters is 

elongation, which can be taken as a substructure 
indicator (although not a definitive one). In most 
cases, clusters with elongated galaxy distributions 
will have substructure. We use Fourier analysis 
(Rhee, van Haarlem and Katgert 1991) to esti- 
mate the elongation of galaxy clusters. In this 
method we assume the azimuthal galaxy distribu- 
tion {N{(j))) resembles the model 



[l + (7Vi/Aro)cos(2<^-2</.o)], (4) 



where A^o is the number of cluster galaxies, 00 is 
the cluster position angle and A^i is the elongation 
amplitude. The statistic for this test is given by 
the elongation strength, defined by 



FE = 



2(52 + (7^)1/2 



\/2Ao 



\/2Ao 



(5) 



where 



11 



S = Y, N{(j))sin{2(f)), C = Y, N{<f))cos{2(P). (6) 

Pinkney et al. (1996) recommend a more strict 
criteria when using FE to detect substructure. 
They advise using FE to reject the circular hy- 
pothesis at the 1% level, with FE being greater 
than 2.5. This is motivated by the fact that elon- 
gation is also a signature of other formation sce- 
narios besides hierarchical mergers. Rather than 
adopt these somewhat arbitrary criteria, we opti- 
mize this test the same way as the others (section 
5.5). The first criterion (1% rejection level) may 
be too rigorous, while the second (selecting only 
clusters with FE > 2.5) makes little difference to 
the final results (Figure 9). 

5.4. The Lee Statistic 

The Lee statistic (Lee 1979) tests for bimodality 
in a given distribution, which is perhaps the sim- 
plest type of substructure. In our case, we test the 
hypothesis that the galaxies in a cluster lie in two 
clumps. Since its introduction in 1979, the test has 
been used in astronomical applications, e.g. Fitch- 
ett & Webster (1987) and Rhee, van Haarlem and 
Katgcrt (1991b). The algorithm works as follows: 
(i) the data arc projected onto a line, which makes 
an angle (j) with a second line, assumed to be of 
constant declination, (ii) the first line is then ro- 
tated by small increments for 0° < (f) < 180°. (in) 
For each orientation all points are projected onto 
the second line, so that each point assumes a new 
coordinate, Xj, along that line, (iv) A search for 
the best partition into two clumps (named "left" 
and "right") is performed. The quantities ai, ar 
and ctt are calculated for the left, right and to- 
tal samples, respectively, for the TV — 1 partitions. 
The a are given by cr = "^{xj — /u)^, with given 
by X^ajj/n. The quantity L is then defined by 

L = maXpartitions ( ; 1 ) • (7) 

\ai+ar J 

L varies with 0, becoming large when the projec- 
tion axis connects two distinct clumps in the orig- 
inal dataset. Following Fitchett (1988) we adopt 
Lrat = Lmax/Lmin = maxL^cj}) / minL{(j)) as this 
test's statistic. As in Pinkney et al. (1996) we use 
6° increments in </>. As pointed out by Pinkney 
et al. (1996), the Lee 2D test is a conservative 



test because it docs not detect elongations as sub- 
structure and also loses sensitivity in the presence 
of more than two subclusters (or if there are two 
clumps, but their sizes are dissimilar). 

5.5. Summary of Substructure Tests 

The codes used for the substructure analysis 
arc those of Pinkney et al. (1996), who evaluated 
the performance of 31 statistical tests (including 
the four two-dimensional tests described above) 
on A''— body simulations of galaxy cluster merg- 
ers. For the two-dimensional tests they recom- 
mend three diagnostics for substructure studies: 
FE, Lee 2D and /?. Based on the simulations, they 
find that the (3 test is the most sensitive, while 
the AST is very insensitive as it requires a dis- 
tinct compact subclustcr to detect substructure. 
The FE test is very sensitive to elongations, but 
they advise setting a threshold for the FE statis- 
tic. The Lee 2D test requires the presence of a 
very distinct subcluster, and is insensitive to more 
than two clumps (or two clumps with very differ- 
ent sizes). Pinkney et al. (1996) point out that Lee 
2D is sensitive only to genuine two-dimensional 
substructure, while FE and (3 are able to find 
substructure in a cluster which is not necessar- 
ily a merger. These conclusions are based on their 
extensive simulations, and we do not have suffi- 
ciently detailed data to independently test their 
assertions. We describe below how we optimize 
the tests and the results when they are applied to 
the DPOSS data. 

Pinkney et al. (1996) conclude with a note say- 
ing that, when galaxies have no redshifts, sub- 
structure could be erroneously detected due to 
the superposition of foreground and background 
groups. While this is certainly true, the large 
amount of two-dimensional data available from the 
most recent sky surveys will suppress this effect. 
The acquisition of deep X-ray data or galaxy ve- 
locities for a large sample of moderately redshift 
galaxy clusters is very expensive, observationally 
speaking. Thus 2D data are also useful for the se- 
lection of subsamples for future surveys (Kriessler 
& Beers 1997). Kolokotronis et al. (2001) used 
an alternate approach when computing the frac- 
tion of optical substructure estimates that are due 
to background contamination, searching for sub- 
structure using 2D optical and X-ray imaging data 
for 22 clusters. They find good agreement between 
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optical and X-ray results, with only ^ 22% of the 
clusters yielding inconsistent results (namely, sub- 
structure seen in the optical, but not in the X- 
ray). They argue that this difference is probably 
due to projection effects in the optical data. How- 
ever, it is important to point out that some recent 
studies have found that projection can also affect 
the morphology inferred from X-ray data (Lopez- 
Cruz et al. 2004; Rines & Diaferio 2006). In our 
case, wc try to minimize this problem by using 
only galaxies selected in a given magnitude range 
(an issue that is not very well addressed in previ- 
ous works). Wc have also optimized our methods 
to minimize the scatter when comparing optical 
and X-ray properties {e.g. Tx—^gais, described 
below). 

5.6. Dependence on Input Data 

The substructure results we show below are 
based on 10190 NoSOCS clusters at 0.07 < 
z < 0.21, which is the largest sample used for 
this type of study to date. The rationale for using 
only clusters at these rcdshifts is detailed below. 
For the optimization of the tests we have used 
the subsample of 445 clusters with a counterpart 
in BAX. Wc first need to quantify the sensitiv- 
ity of the cited tests to the choices of centroid, 
redshift, magnitude range, maximum radius and 
significance level by evaluating substructure for 
different ranges of these parameters. In Figure 
6 we show the effects of the adopted center and 
redshift. We test four centers, the original cluster 
position (O), the BAX coordinate (B), the lumi- 
nosity weighted center (L) and the position of the 
maximum density within the cluster (P). The six 
panels in this figure are analogous to those in Fig- 
ure 1, now showing the residuals in the /3 test. 
Additionally, we show in the upper right of each 
panel the value of IQR for the residuals of the 
significance of the /3 test (not plotted; these are 
labeled Qcr2)- The largest differences are found 
when comparing the original cluster position to 
the BAX coordinate (or the density peak posi- 
tion) and when adopting different redshifts (BAX 
or DPOSS). The centroid choice affects the general 
galaxy distribution, which may not be well indi- 
cated by the X-ray center or maximum density 
peak (see the comments regarding Figure 1). The 
redshift affects both the sampling of the cluster 
luminosity function and the apparent radius used 



for selecting galaxies for the substructure tests. 

Based on this figure and on Table 1 we evaluate 
substructure (for the entire sample of 10190 clus- 
ters) using the original cluster position, instead of 
the luminosity weighted center or the density peak 
coordinate. For consistency, when using substruc- 
ture to exclude clusters in comparisons of optical 
and X-ray properties (§6) we also adopt the orig- 
inal cluster position and the BAX redshifts (see 
the richness discussion). As shown in Table 1 this 
gives the minimum scatter for the optical— X-ray 
relations based on the substructure results. 

In Figure 7 we show the dependence of the sub- 
structure results on the magnitude range (lower 
panel) and maximum radius (middle panel) used 
when selecting projected cluster galaxies. Due 
to the large computational cost of running these 
tests for the entire sample of NoSOCS clusters, we 
show the dependence on the magnitude range and 
maximum radius only for the STibsct common to 
BAX. Obviously, the ability to detect substruc- 
ture decreases when considering only the bright 
central galaxies. Conversely, at fainter flux limits 
and larger radii, the 2D data presents problems 
due to background and foreground contamination. 
Wc first attempted to estimate substructure us- 
ing all galaxies within the survey magnitude lim- 
its (15. 0<r<19. 5), but this leads to misleading 
results due to different sampling of cluster galax- 
ies for clusters at different redshifts. We thus in- 
vestigated the dependence of the number of clus- 
ters with substructure on the magnitude range 
used. We tested 6 ranges as indicated in the lower 
panel of Figure 7, using only 100 (out of 638) clus- 
ters that span the full magnitude range sampled 
(m* — 2 < rrir < m*+2, corresponding to redshifts 
of 0.11 < z< 0.14). 

We conclude that the high percentage of clus- 
ters with substructure when using faint galaxies 
{m*+2) may be an artifact due to background con- 
tamination. Wc do not require the brightest galax- 
ies (to* — 2) because that would further reduce 
our sample, eliminating lower redshift clusters. 
As shown in Figure 7, the results for the ranges 
TO* — 2 < rrir < TO* -|- 1 or to* — 1 < rrir < to* -|- 1 
are similar, having little effect on our conclusions. 
For all substructure analyses we therefore use the 
latter range, which is fully sampled by our data 
at redshifts 0.07 < z< 0.21, where we find 445 of 
the 638 common BAX and NoSOCS clusters (430 
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Fig. 6. — The (3 test residual for different cluster centroids and adopted redshifts. The comparison shown in 
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we show in the upper right of each panel the values of IQR (named Qo-2) for the residuals of the significance 
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of the 618 with estimates of Lx and NgaZs) and 
53 of 101 with estimates of Tx- From the 16546 
clusters in NoSOCS there are 10190 in this redshift 
range. 

From the middle panel of Figure 7 wo see that 
out to R = 1.00 Mpc the substructure tests 
behave similarly, except for AST, which is recog- 
nized by Pinkney et al. (1996) as being very insen- 
sitive. In the upper panel of Figure 7 we show the 
dependence of the number of clusters with sub- 
structure on the significance level (S.L.) required 
to reject the null hypothesis. Obviously, as we re- 
lax the threshold to consider a substructure esti- 



mate meaningful, the number of selected clusters 
increases. Our choice of significance level is ex- 
plained below. 

Any optical study relying on purely projected 
galaxy positions is affected by superpositions of 
loose groups. One could assess projection effects 
in substructure estimates by using X-ray data 
(Kolokotronis et al. 2001) or including velocity dis- 
persion information. Given the large number of 
clusters in our sample, this is not practical. Since 
our goal is to compare optical and X-ray cluster 
properties (§6) we optimize the parameters (mag- 
nitude range, radius and S.L.) for the substructure 
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tests by minimizing the scatter in these scahng re- 
lations. This could, in principle, lead to similar 
results as those obtained by Kolokotronis et al. 
(2001). This strategy is also complementary to 
the analysis in Figure 7, which gives no objective 
criteria to choose the optimal parameters. 

This approach is similar to the one show in Fig- 
ure 2. Before, we wanted to find the optimal radius 
to compute richness (the one that minimizes the 
scatter of the scaling relations). Now we want to 
find the optimal radius, magnitude range and S.L. 
to estimate substructure. Note that the optimal 
radius for substructure measurement need not be 
equal to the one used for richness, as we are mea- 
suring different cluster properties. 

In the upper panels of Figure 8 we examine 
how the orthogonal scatter of three correlations 
of cluster parameters (J-ix—^gais, T^x—^gais and 
Tjf— Ljf) can be minimized by excluding clusters 
with signs of substructure. In the left panel, we 
show the dependence of the scatter on the radius 
used for selecting cluster galaxies, while on the 
right we show the variation with the significance 
level employed. The results are shown only for 
the (3 test. As each test has a different depen- 
dence on these parameters, their optimal values 
will also differ, implying that two tests could yield 
the same percentage of clusters with substructure 
through the use of very different apertures and es- 
pecially significance levels. After many trials for 
all four tests we find that the (3 test is most appro- 
priate (for our data) for detecting substructure. It 
gives results that are comparable to the literature 
and is considered the most sensitive among the 
four tests (Pinkney et al. 1996). Based on Figure 
8, we use a radius of 1.50 h^^ Mpc and a sig- 
nificance level of 5%. A similar analysis for the 
magnitude range suggests that our original choice 
(m* — 1 < nir < m* -|- 1) is appropriate. 

It is interesting that only the Tx—^gals rela- 
tion is sensitive to the choices of radius and S.L. 
As noted in Section 6, the scatter in the Lx—^gais 
relation is much higher than for Tx—^gals- The 
exclusion of clusters based on substructure is not 
sufficient to improve the l^x—^gais relation. In 
the lower panels of Figure 8 we show similar plots, 
but for the four substructure tests and only for the 
Tx—^gais relation, using our final settings of R = 
1.50 Mpc and S.L. = 5%. From these two 
panels it is clear that only the /3 test is sufficiently 
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tical tests as a function of the magnitude range 
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panel); the significance level required for each test 
(upper panel). 



sensitive. 

The distribution of substructure results, for the 
four statistical tests, is shown in Figure 9. The 
solid lines represent the distribution of all clusters, 
while the thick dotted lines show the subset with 
substructure. A cluster is considered to have sub- 
structure if the null hypothesis is rejected at the 
5% significance level. This means that at most 5% 
of the Monte Carlo datasets yield a statistic more 
extreme than the real cluster. 

We have also investigated the possibility that 
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the substructure tests depend on the cluster con- 
trast, richness and redshift, with the results shown 
in Figure 10. As mentioned in §4.1, the con- 
trast is defined as the ratio between the number 
of galaxies in the cluster region (within the same 
magnitude range sampled for the richness calcula- 
tion) and the background error Qa^kg ■ There is no 
clear sign that substructure detection is contrast- 
dependent. The variation in the fraction of clus- 
ters showing substructure with the cluster con- 
trast is shown in lower panel, the redshift vari- 



ation in the middle panel, and the richness depen- 
dence in the upper panel. We can not verify a 
clear trend with contrast, but see a mild increase 
with richness and a very small decrease with red- 
shift. Since this analysis is restricted to clusters 
at 0.07 z ^ 0.21, it is difficult to investigate the 
evolution of substructure with look-back time. We 
postpone the study of the evolution of substruc- 
ture with redshift to a future work, where we will 
use deeper photometry (reaching higher redshifts) 
combined with X-ray imaging. 
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Fig. 9. — Distribution of the substructure results for the four statistical tests: (a) the angular separation 
test (AST); (6) the /3 test (to improve the visualization the (3 values are multiplied by 1000); (c) Fourier 
Elongation (FE) and {<£) Lee 2D results. The solid lines represent the distribution of all clusters, while the 
thick dotted lines show the results for the subset with substructure based on the Monte Carlo simulation, 
rejecting the null hypothesis at the 5% significance level. Note that AST is the only test for which lower 
values represent a more extreme statistic. For the FE test wc show the distribution of all clusters considered 
to have substructure at the 5% significance level, without further excluding clusters with FE > 2.5. As we 
can see the adoption of this additional criterion - as suggested by Pinkney et al. (1996) - excludes only a 
small additional number of clusters. 



The final substructure results using the 10190 
clusters at 0.07 ^ z ^ 0.21, for each test are as 
follows: AST finds that 13% clusters have sub- 
structure; the /? test selects 35%; the FE test finds 
21%; and the Lee 2D test finds 45%, all at the 
5% significance level. As we can see, AST selects 
by far the fewest clusters with substructure, as 
expected due to its lack of sensitivity. The rates 
found by the /? and Lee 2D tests are in good agree- 



ment with the results of Kolokotronis et al. (2001), 
who found at least ^ 45% of their clusters with 
strong indications of substructure. However, their 
sample was much smaller (only 22 clusters), and 
they employed other tests, using optical and X- 
ray data to confirm the substructure estimates. 
Other literature results indicate that the percent- 
age of clusters showing substructure varies from 
30% to 80%. This discrepancy is generally be- 
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lieved to be associated to the method employed 
for each analysis (Kolokotronis et al. 2001). We 
stress that the variation in these results is likely 
also due to the cluster centroid, area, magnitude 
range and significance level adopted when apply- 
ing a given test. For instance, had we used 1% 
or 10% as our significance level for the /? test, 
we would find 19% and 46% of clusters with sub- 
structure, respectively. Thus, not only it is funda- 
mentally important to compare similar datasets, 
the tests employed must be consistent. The same 
clearly holds true for comparisons with cosmolog- 
ical simulations. 

5.7. Correlation Between Substructure 
and EUipticity 

EUipticity is generally considered to be a sign of 
substructure, a result that has been found by pre- 
vious authors (Jones & Forman 1999; Kolokotro- 
nis ct al. 2001). We thus decided to investigate 
if the substructure tests we employ in this work 
are correlated to ellipticity. For this purpose we 
used only the 445 common BAX-NoSOCS clus- 
ters at 0.07 ^ z < 0.21. In order to compare our 
results to the ones from Kolokotronis et al. (2001), 
we have considered siibstructurc mcasiires within 
0.75 Mpc, an aperture close to the maximum 
radius employed by these authors (0.60 Mpc). 
We have also binned the substructure and ellip- 
ticity results into substructure bins. For the four 
tests used in this work we fit a linear correlation of 
the type e = A -I- Bsubjest , where e is the ellipticity 
value, subtest is the result for a given substructure 
test, and A and B are the intercept and slope ob- 
tained by an orthogonal regression method (see 
§6.2). The relations obtained for all the four sub- 
structure tests are shown on Figure 11 and sum- 
marized in Table 2. We note that two of the meth- 
ods (FE and Lee2D) show a clear dependence be- 
tween ellipticity and substructure. These results 
are in agreement to the ones found by Kolokotro- 
nis et al. (2001), although we find slopes shallower 
than what they found. However, it is important to 
keep in mind that they employed a different sub- 
structure test and selected galaxies from a smaller 
aperture and difii'erent magnitude range. 

As can be seen from Figure 9, the /? test can 
take on negative values. We found that these are 
anti-correlated to e, while the positive values are 
correlated. We have thus considered only the pos- 
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Fig. 10. — Substructure results for the /3 test as a 
function of cluster contrast (lower panel), redshift 
(middle panel) and richness (upper panel). 



itive values of /3 for the results shown in Table 2 
and Figure 11. From inspection of this table we 
see that AST is the only test that yields a slope 
consistent to zero. This result is consistent with 
the aforementioned lower sensitivity of this test 
(Figure 7). 

6. Optical vs. X-ray Cluster Properties 

In this section we investigate the correlation be- 
tween richness {^gals) and X-ray properties (Lx 
and Tx) of clusters common to DPOSS and BAX. 
From the 638 common clusters, 620 have Lx (of 
which 2 have no richness estimate), while Tx is 
available for 101 clusters. As mentioned before, 
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Fig. 11. — Correlation between ellipticity and substructure for the four statistical tests: (a) the angular 
separation test (AST); (6) the /3 test; (c) Fourier Elongation (FE) and {d) Lee 2D results. 



for analyses involving substructure estimates we 
restrict the study to 0.07 ^ z ^ 0.21, reducing the 
original samples from 618 to 430 clusters (with Lx 
and richness determined) and from 101 to 53 sys- 
tems (with Tx)- Clusters with Ngais < 10 are also 
excluded when establishing these correlations (see 
the discussion in §4.1 and Figure 3) and a 3— ct 
clipping is also applied to exclude outliers. We 
also discuss how the scatter in the established rela- 
tions is affected by the radius adopted for calculat- 
ing richness, the redshift difference (Az) between 
BAX and DPOSS, positional offsets, the contrast 
of the cluster, and the presence of substructure. 



6.1. Dependence on Input Data 

We reiterate that the orthogonal scatter (cr) and 
the linear correlation coefficient p of the Lx—^gais 
and Tjf — Ngais relations are extremely sensitive to 
the aperture used for calculating richness (see Fig- 
ure 2). These results are in good agreement with 
those of Popesso et al. (2004), who verified the de- 
pendence of the scatter in the Lx—l^opt relation on 
the aperture used. When using very small aper- 
tures to estimate N^ais the galaxy number count 
becomes very uncertain due to centroid errors. 
For large radii the background contribution grows 
rapidly, contaminating the richness estimates. As 
explained in §4.1 we use a 0.50 Mpc radius for 
estimating l^gais for comparisons to X-ray proper- 



19 



Table 2 

Correlation between ellipticity and substructure. 



Subst. Test 


Linear Corr. Coef. 


Intercept 


Slope 


P 


0.92 


0.358 ± 0.017 


0.092 ± 0.052 


Lee2D 


0.94 


0.023 ± 0.037 


0.149 ± 0.017 


FE 


0.98 


0.105 ± 0.014 


0.126 ± 0.007 


AST 


-0.59 


0.415 ± 0.025 


-0.004 ± 0.001 
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Fig. 12. — Variation of the scatter of the Lx—^gais relation as a function of the maximum absolute redshift 
difference between the NoSOCS and BAX values (left panel) and the positional offset (right panel). 



ties. 

The effects of Az and positional offset on the 
scatter in the l^x^^gais relation arc shown in the 
left and right panels of Figure 12, respectively. 
As before, we adopt the BAX redshift to calcu- 
late Ngais, since the richness is strongly rcdshift- 
dependent (through the apparent radius and the 
magnitude range). This should minimize the scat- 
ter in the hx^^gais relation, and we expect this 
scatter to show no dependence on Az if the red- 
shift errors are random. The variation shown in 
the left panel of Figure 12 is small (~ 0.020) but 



real, and is probably due to the fact that the BAX 
redshifts are not homogeneous measures. They are 
automatically retrieved from NED, and are most 
often spectroscopic, but some are photometric es- 
timates. Even the spectroscopic redshifts may be 
based on only one or two galaxy spectra. Thus, 
we conclude that even after adopting the BAX 
redshift, we could further reduce the scatter in 
the Lx^'^gais relation by restricting the sample 
to clusters with small values of Az (at the cost 
of reducing the sample size). The improvement in 
the least-square fit of the Lx—^gais relation for 
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clusters with small positional offset is expected. 
The dependence of (Jorthogonai is of the same order 
as that shown in the left panel. 

We also investigate the dependence of the 
Lx— Nga(s with the minimum and maximum red- 
shifts used to select clusters for the fit. The results 
are optimal for clusters at 0.15 < z < 0.30, with 
strong sensitivity to the lower redshift limit. This 
is because most of the low richness systems are 
at low redshift, and these poor clusters have the 
largest richness errors (Figure 3). Based on Fig- 
ures 2 and 12 we use a radius of 0.50 Mpc 
for calculating richness, but do not apply any 
cut with redshift or offset. All the dependences 
of (Jorthogonai ou the rcdshift cutoff, A^; and the 
positional offset are much less pronounced than 
the dependence on aperture. Thus, we do not re- 
strict the sample based on these weaker effects, 
which would significantly reduce the number of 
clusters and make it difficult to quantify the most 
significant effects. The impact of contrast and 
substructure are investigated further below. 

6.2. Dependence on Fitting Method 

Before examining the connection between op- 
tical and X-ray parameters wc investigate which 
linear regression solution is more robust, using the 
packages slopes and bees-regress from Isobe et al. 
(1990) and Akritas and Bcrshady (1996), respec- 
tively. The first performs unweighted fits, while 
the second considers errors as well as intrinsic scat- 
ter. Because BAX provides no error estimate for 
Lx, we perform an unweighted fit (using slopes) 
for the Lx—^gais relation, while bees-regress is 
used for Tx^'^gais- Wc use two commonly used 
solutions for each package. The first is the ordi- 
nary least-square (OLS) bisector solution, the line 
which bisects the OLS(Y|X) (minimizing the resid- 
uals in Y) and OLS(X|Y) (minimizing the residu- 
als in X) lines. The second solution is the orthog- 
onal regression line, which minimizes the perpen- 
dicular distances (Isobe et al. 1990; Feigelson et 
al. 1992; Akritas and Bershady 1996). We stress 
the comment of Isobe et al. (1990) concerning lin- 
ear regression methods that the two methods used 
here (and others available in the literature) find 
regression coefficients that are theoretically differ- 
ent from each other. In other words, they do not 
represent different estimates of the same quantity. 
Even if the entire population were sampled, the 



measured slopes would differ (only in special cases 
is a single relation obtained, such as when p= 1). 

In order to evaluate which of these methods is 
best suited for our purposes, we generated artifi- 
cial data sets with known scatter and slopes. The 
relations resulting from the optical - X-ray com- 
parison could be very steep and noisy, especially 
those involving hx- Hence, we created data sets 
of 1000 points for different angles (ranging from 1 
to 89 degrees) and for 20 different scatters (from 
0.01 to 0.40). We then nm slopes and bees-regress 
for the bisector and orthogonal solutions, comput- 
ing the difference between the input and output 
slopes (lAgjopel)- This procedure is repeated 100 
times. The results are shown in Figure 13, where 
the error bars for each point represent the stan- 
dard deviation for the 100 events. In the upper 
panels we show the variation of |As(ope| with angle 
(at a fixed scatter of 0.20), while the dependence 
with scatter is shown in the lower panels (for a 
fixed angle of 65°). The tests using unweighted 
and weighted fits are shown in the left and right 
panels, respectively. The bisector results are plot- 
ted as solid lines while the orthogonal solutions 
are the dotted lines. As with the real data, a 3— a 
clipping is applied to exclude outliers. The orthog- 
onal minimization is clearly superior for recover- 
ing the input slopes, especially when there is large 
scatter. If we do not eliminate outliers, the results 
shown in the upper panels diverge for large angles, 
but the bisector solution still shows larger devia- 
tions from the original slopes. Given these results, 
we use the orthogonal regression for all fits. We 
note that some results in the literature (Marke- 
vitch 1998) have used the bisector, which may be 
misleading. 

6.3. The X-ray and Optical Relations 

6.3.1. X-ray Luminosity and Temperature 

We begin our analysis by examining the Lx— Tx 
relation (Figure 14) derived from the BAX data, 
to verify the reliability of this data source for our 
purposes. We use 120 clusters with available tem- 
peratures within the NoSOCS region. As the solid 
thin line we show the best fit using all 120 clus- 
ters (two are excluded as outliers). As the thin 
dotted line we show the results obtained using 
only the 100 clusters with optical contrast > 10.0. 
As discussed below, the scatters and slopes of 
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Fig. 13. — Variation of |As;ope|, the difference between tlie input slope and the one recovered by the bisector 
and orthogonal solutions. In the upper panels we show the variation with input angle (at a fixed scatter of 
0.20), while in the lower panels wc show the variation with scatter (at a fixed angle of 65°). The left panels 
show the results for the unweighted fits and the right panels for the weighted fits. The bisector results are 
given by the solid lines, while the orthogonal ones are shown by the dotted lines. 



the scaling relations depend on the contrast cut 
applied to the cluster sample. Interestingly, the 
optical contrast appears to provide a useful test 
for selecting clusters to affect the Lx—Tx rela- 
tion. Literature results are shown by thick lines. 
The solid and dotted lines show the results from 
Markcvitch (1998), with the dotted line for Lx.boi- 
We note that Markevitch (1998) employed a bisec- 
tor method for fitting this relation. If we also use 
the bisector slope, our results arc considerably less 
steep (Ljf oc Tx^'"^^), showing better agreement 
with Markevitch (1998). However, as discussed 
above, the bisector is not the optimal method for 



fitting steep slopes especially in the presence of 
large scatter. The results from David et al. (1993) 
are plotted as the dotted thick line. 

6.3.2. X-ray Luminosity and Richness: Struc- 
ture Segregation 

The Lx—^gais relation is shown in Figure 15, 

along with the effects of the contrast and the pres- 
ence of substructure in the derived relation. The 
comparison of all clusters at 0.07 <i;< 0.21 is ex- 
hibited in the lower left panel, while the high con- 
trast systems (contrast > 12.0) are shown in the 
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Fig. 14. — The Lx - Tx relation obtained using all 120 BAX clusters (with two excluded as outliers) in the 
NoSOCS region is shown as a thin solid line. The sample is then restricted to high optical contrast systems 
(C > 10.0) with the result indicated by the thin dotted line. Literature results are indicated by thick lines. 
The relations obtained by Markevitch (1998) are indicated by the solid and dotted lines (the last one uses 
bolometric luminosities). The dashed Hne indicates the results from David et al. (1993). 



lower right. The substructure-free clusters are ex- 
hibited in the upper left panel while those with 
substructure in the upper right. The fit shown by 
a solid line in all panels is that obtained for the 
substructure-free clusters. In panels (b), (c) and 
(d) we also show the 3(t boundaries of this fit (as 
dotted lines) and the fit obtained with the data 
shown in each panel (all clusters, high contrast 
systems and clusters with substructure) with the 
dashed line. In all panels we indicate the number 
of clusters used in the fits (after excluding outliers 
at a 3— (T level), the linear correlation coefficient 
(p), the intercept (A), slope (B) and the orthogo- 



nal scatter of the relation obtained for that sam- 
ple. 

From this figure we see that the scatter changes 
only slightly when selecting only substructure-free 
clusters, while it is reduced by ~ 12% when se- 
lecting high contrast systems. The combination of 
the two criteria (substructure and contrast) does 
not improve the scatter compared to applying only 
the contrast cut. Interestingly, the slope is sig- 
nificantly higher for the high contrast clusters. 
Popesso et al. (2004) also found improvement in 
the relations when excluding low contrast systems. 
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Wc conclude that the scatter in the hx - ^gais 
relation is mainly due to the difSculty of obtain- 
ing accurate measures of these parameters for low 
luminosity systems. If there is any difference for 
clusters with or without substructure, it is washed 
out in these relations. Thus, when correlating Lx 
to the optical parameters wc do not detect any 
significant structure segregation. 

6.3.3. X-ray Temperature and Richness: Struc- 
ture Segregation 

The Tx— Ngais relation is plotted in Figure 16 
(with 3— a outliers shown as open symbols), in a 
manner analogous to the previous Figure. In the 
lower left panel of Figure 16 we show the results 
for the 53 clusters at 0.07 < ^ < 0.21 with Tx and 
richness measiircd. As before, in the lower right 
panel we show the results for the high contrast sys- 
tems (contrast > 12.0). In the upper left the sam- 
ple includes only clusters found to have no s\ib- 
structure by the /3 test at the 5% S.L, while in the 
upper right we show the complementary clusters. 
A few things stand out from this figure. First, 
Tx is clearly better correlated with the optical 
parameters than Lx- Second, substructure has 
an effect on the scatter of the Tx—'^gais relation, 
which was not seen when using hx ■ When exclud- 
ing clusters with substructure the scatter of the 
Tjsf— Nga(s relation is reduced by ^ 63%. Third, 
there is no improvement when applying the con- 
trast cut. This is mainly due to the fact that Tx is 
typically available only for the most massive sys- 
tems, so the analysis is restricted to a small range 
of temperature (and also contrast). Of the 430 
clusters with Lx measured, 177 (^ 41%) clusters 
are low contrast systems, while they are only 14 
(~ 26%) of the 53 with T^ . Finally, when com- 
bining the two criteria (substructure and contrast) 
the scatter is not very different than if we employ 
only the substructure cut, implying that the dom- 
inant factor in the scatter of this relation is sub- 
structure. 

The small number of clusters with available 
temperatures is certainly a problem when estab- 
lishing these relations and applying the above cuts. 
However, Smith et al. (2005) have also recently 
foimd that the scatter in the scaling relations be- 
tween cluster mass. X-ray luminosity and temper- 
ature is dominated by unrelaxed clusters. Their 
sample is smaller than ours (and they are also lim- 



ited to massive systems), but serves to corroborate 
the current results. An opposite result (no struc- 
ture segregation in the scaling relations) is found 
by O'Hara et al. (2006), using both observational 
data and simulations. They compare the devia- 
tions of clusters from the best fit relations (ob- 
tained for all systems and not for the separate sets 
with or without substructure) for different scaling 
relations and find no evidence for higher scatter 
for clusters with more substructure. The substruc- 
ture tests they employ are different than the ones 
used here. Their method for dividing the sam- 
ple into clusters with or without substructure is 
also unclear (there is no cut at a given significance 
level). However, the most important differences 
between our methodology and theirs are the fact 
that they investigate substructure with an aper- 
ture that scales with the cluster mass (-R500 and 
-R2500), while we use a fixed metric (1.50 Mpc 
radius), and that they rely on X-ray data, while 
we use optical data. Clusters with more substruc- 
ture arc probably dynamically younger and would 
thus be expected to deviate from relations estab- 
lished for well relaxed systems. In future work, we 
will investigate siibstructurc using these alterna- 
tive approaches (O'Hara et al. 2006) to see if our 
current results persist. 

It comes as no surprise that the X-ray tem- 
perature correlates better with optical properties 
than the X-ray luminosity. Tx has been shown to 
provide a much tighter correlation to cluster mass 
than X-ray luminosity (Voit 2005), as Tx is closely 
related to the cluster's potential well. What we 
find here is that the connection between richness 
(which is a mass tracer) and Tx has a smaller 
scatter than if we used Lx. Similar results are 
foimd by Yce and EUingson (2003) and Popcsso et 
al. (2005) who find tighter relations between their 
richness estimates (Bgc or Lopt) and Tx than with 
Lx ■ Wc also note that samples constructed purely 
at one wavelength may not include the same types 
of objects (Lubin et al. 2004). 

6.3.4. X-ray Luminosity and Richness: Sample 
Dependence 

We now compare our findings to literature re- 
sults. We show in Figure 17 the connection be- 
tween Lx and richness or Lopt. Our results are 
shown by thin lines, while the literature ones are 
thick lines. When performing this comparison we 
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Fig. 15. — The Lx—^gais relation is shown for: (a) aU clusters at 0.07 ^ z < 0.21; (6) clusters with contrast 
> 12.0; (c) clusters without substructure (at the significance level of 5%) based on the /3 test; (d) clusters 
with substructure according to the (3 test. The fit shown in all panels (as a solid line) is the orthogonal 
solution obtained for the sample in the upper left panel. For three panels (a, b and d) the dotted lines show 
the 3a boundaries for this fit. For these 3 panels (a, b and d) the fit obtained with the data shown on each 
panel (all clusters, high contrast systems and clusters with substructure) is shown by the dashed line. 



investigate a few effects: (i) how the slope of the 
optical vs X-ray relation changes with cluster con- 
trast; (ii) the impact of the ccntroid (optical or 
X-ray) on the slope; (iii) the slope difference when 
using richness or Lopt as the optical parameter; 
(iv) the difference in the results obtained by the 
bisector and orthogonal fits. All of the relations 
based on our data shown here include all clusters 
(we do not exclude clusters with NgaU < 10.0, 



the only exception being the relation for clusters 
with no substructure) and are obtained using an 
orthogonal regression. We also stress that the rela- 
tions based on the X-ray center (BAX) are derived 
from all BAX clusters within the area covered by 
NoSOCS (whether or not they have a counterpart 
in DPOSS). The relations obtained with the op- 
tical center include only clusters common to the 
optical and X-ray catalogs. Hence, the former re- 
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Fig. 16. — The Tx-^gais relation for all clusters at 0.07 ^ z < 0.21. The full sample is shown in the lower 
left panel; high contrast systems are exhibited in the lower right panel; clusters considered to be substructure 
free arc exhibited in the upper left panel, where two outliers are shown with an open symbol; clusters with 
substructure are shown in the upper right panel. 



lations are based on larger samples of clusters. 

Figure 17 shows the Lx—^gais relations for a 
variety of samples. All of the results from our 
sample are plotted as thin lines, including: 

• Solid line: all BAX clusters in the NoSOCS 
region. Of the 792 BAX clusters, 744 have 

both richness and Lx available. Ten clus- 
ters are eliminated as outliers after a 3— a 
clipping is applied. 



• Dotted line: The previous sample trimmed 
to clusters with optical contrast > 5.0. 

• Short dashed line: The sample is further re- 
duced to clusters with contrast > 10.0. For 
the above three relations we use the X-ray 
centroid when estimating 'Ngais ■ 

• Long dashed line: As the solid line above, 

but using only the common clusters between 
NoSOCS and BAX (with available Ljc and 
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Ngafs), and the optical centroid for the rich- 
ness calculation. 

• Dashed— dotted line: The sample is then cut 
to consider only substructure free systems, 
using the optical centroid (as shown in Fig- 
ure 15). 

• Long dashed— dotted line: The Lx—l^opt re- 
lation, using the X-ray centroids. 

We also show a variety of results from the litera- 
ture as thick lines: 

• Solid line: The hx - ^opt,i relation obtained 
by Popesso et al. (2004) for their full sample. 
We note that they too use the X-ray centers. 

• Dotted line: The relation from Popesso et 
al. (2004) when groups are excluded. 

• Short dashed line: The correlation between 
Lx (bolometric) and optical richness (given 
by Kci) from Donahue et al. (2001). 

• Long dashed line: The correlation between 
hx and richness (Bgc) obtained by Yee and 
Ellingson (2003). 

We note that Yee and Ellingson (2003) use a bi- 
sector solution, which is generally less steep. If we 
considered a bisector solution in our analysis the 
powers in the relations we obtained (thin lines) 
would range from 1.79 to 2.19, which is in good 
agreement with their results. 

From this figure we can sec that the slope in- 
creases rapidly as we consider higher contrast sys- 
tems. This trend was already evident in Figure 
15, where we also sec the scatter is greatly reduced 
when the sample is restricted. However, this be- 
havior does not persist for much larger values of 
the contrast cut. In Figure 17 we show the results 
for the full sample and with cuts at C > 5.0 and 
10.0. If we increase the cut to C > 15.0 or 20.0 
the slope variation ceases, with the slope actually 
becoming lower. That is due to sampling effects as 
the number of clusters is drastically reduced. It is 
important to mention that contrast is directly pro- 
portional to richness. Thus, if richness cuts were 
applied - instead of contrast cuts - we would find 
similar trends. An important issue is that when 
using low luminosity or richness systems the data 



is likely to be censored. Donahue et al. (2001) ex- 
tensively discuss this problem and show how the 
correlation between Lx and richness (Ac;) changes 
when the fit is derived in conjunction with the A,,; 
function and the X-ray luminosity function [i.e., 
taking into account sampling effects at low flux 
regimes). 

Based on Figure 17, a very interesting conclu- 
sion is that the relations are less steep when us- 
ing the X-ray centroid (instead of optical) to cal- 
culate richness (or Lopt). When using all clus- 
ters and considering the X-ray center, the slope 
is 3.73, increasing to 4.29 and 5.60 when trim- 
ming the sample at C > 5.0 or 10.0, respectively. 
The slopes measured when considering the opti- 
cal centroid for these three cases (using the com- 
mon BAX and NoSOCS clusters) are 4.54, 5.09 
and 6.40. This effect is associated to the contrast 
cuts, as the galaxy number density is generally 
reduced when we consider the X-ray centroid, in 
which case we expect the slopes to be smaller. We 
also find that the relation using optical luminos- 
ity is shallower than the one with richness. When 
comparing our results to the literature we note 
that Popesso et al. (2004) found the same effect 
regarding the contrast cut (by excluding groups 
from their main sample) . The relations we obtain 
are also in good agreement to their results, as well 
as those from Donahue et al. (2001) and Yee and 
Ellingson (2003). In the latter case there is con- 
cordance if we consider only the bisector solution. 

We define ^gais differently from the richness 
parameters in the literature shown here (Ac; and 
Bgc). Thus, the most direct comparison of our 
results to the literature is made when consider- 
ing the relations involving Lopt. The results from 
Popesso et al. (2004) show less power than ours 
(2.13 in comparison to 2.54). However, their sam- 
ple may be composed of a larger number of low 
luminosity systems (groups) than ours. When 
they exclude groups, the power increase to 2.78, 
in much better agreement with our findings. Even 
considering their full sample our results agree at 
the 2(7 level, while the agreement is at the la level 
for their high contrast sample. The use of the 
optical center (instead of the X-ray centroid) for 
computing Lopj yields steeper relations. When us- 
ing all the common BAX-NoSOCS clusters we find 
a slope of 3.01 for the Lx-hopt relation, which dis- 
agrees with the full sample of Popesso et al. (2004) 
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Fig. 17. — The correlation between Lx and richness or I^opt- The thin lines show the results obtained in this 
work, while the thick lines exhibit results from the literature. This work: (i) solid line: Lx—^gaU for all 
BAX clusters in the NoSOCS region; (ii) dotted line: the previous sample trimmed at contrast > 5.0; (iii) 
short dashed line: the sample is further cut at contrast > 10.0. For relations i-iii we considered the X-ray 
centroid when estimating Nga/s; (iv) long dashed line: as (i) but considering all common clusters between 
NoSOCS and BAX and using the optical centroid for richness calculation; (v) short dashed— dotted: 
sample (iv) cut to consider only substructure free systems; (vi) long dashed— dotted line: The connection 
between Lx and Lopt (BAX coordinates are used for computing Lgpt) for all BAX systems in the NoSOCS 
region. Literature results: (i) solid line: The Lx ~ Lopt.i relation obtained by (Popesso et al. 2004) for 
their full sample; (ii) dotted line: The relation obtained by the same authors when excluding groups; 
(iii) short dashed line: Lx (bolometric) vs. optical richness (Ac;) from Donahue et al. (2001); (iv) long 
dashed line: the connection between Lx and richness (Bgc) as obtained by Yee and EUingson (2003). 



at the 3(7 level, but is in concordance - at the 2a 
level - to their results excluding groups. 

Another effect that merits attention is the fit- 
ting method employed (see §6.3). Literature re- 
sults are often contradictory, and one of the causes 
that is rarely discussed is the linear regression em- 
ployed. We have shown (Figure 14) that our re- 



sults differ from Markevitch (1998), mostly be- 
cause they adopt a bisector solution when fit- 
ting the Lx-Tx relation. On the other hand our 
findings are consistent with those of David et al. 
(1993). Similar differences are shown when com- 
paring X-ray luminosity and richness. The smaller 
power obtained by Yee and EUingson (2003) in 
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comparison to ours is mainly associated with the 
fitting method (bisector in their case). It is nec- 
essary to consider all of the possible effects (cen- 
troid adopted, flux limit and fitting method, at a 
minimum) when attempting to calibrate a mass- 
observable relation. Systematic errors in the mea- 
surement of cluster mass arc exponentially ampli- 
fied by the steepness of the cluster mass function. 
In future work we will link richness to mass and 
then quantify the effects of different solutions on 
the mass function estimates. 

6.3.5. X-ray Temperature and Richness: Sample 

Dependence 

In Figure 18 we show the connection between 
Tx and richness or Lopt- As before, our findings 
arc shown by thin lines, while thick lines show 
the literature results. We have also not excluded 
clusters with Ngais < 10.0 (except for the rela- 
tion of clusters without substructure), and the 
fits are obtained through orthogonal regression. 
The results shown here include only clusters at 
0.07 < z < 0.21. The thin lines (our work) are 
analogous to those in the previous figure, but for 
the relations involving Tx- The literature results 
(thick lines) arc the Tx—^opt.i relation obtained 
by Popesso et al. (2004) (solid line), the Tjf-Bgc 
relation from Yee and Ellingson (2003) (dotted 
line), the Tx— Ljv/tj relation from O'Hara et al. 
(2006) (short dashed line) and the Tx—^opt,r rela- 
tion obtained by Plionis et al. (2005) (long dashed 
line). We do not find the same trend with contrast 
as for Ljf— Nga(s. This may be due to the smaller 
fraction of clusters excluded based on contrast for 
the Tjjf— Ngajs relation (see above), such that this 
relation is already restricted to high contrast sys- 
tems. There is still a small difference between the 
results with different ccntroids, and "Lopt again de- 
fines a shallower relation than Ngais ■ Our findings 
also show good agreement with the literature re- 
sults. 

A few comments on our results are still nec- 
essary. If it is assumed that mass traces optical 
light and the gas is in hydrostatic equilibrium, wc 
expect the mass- luminosity {M/Lopt) ratio to be 
constant and the gas temperature to be propor- 
tional to mass (T oc M2/3). The X -ray luminos- 
ity is also related to temperature {Lx oc T"), al- 
though the precise value of a is subject to some 
debate. The theoretical expectation is a = 2, but 



its observed estimate is closer to 3 (this work and 
David et al. (1993)). Considering these relations 
and assuming a = 3, for instance, we expect the 
X-ray and optical luminosities to be related as 
{Lx oc ^opt) and the relation to temperature to 
be Tx oc L^pf. If we consider that 'Ngais is directly 
proportional to Lopt our results point to steeper 
relations. Lx oc L^^^ and Tx oc L^pf, although 
Tx—Lopt becomes less steep if we exclude clusters 
with substructure. As pointed out by Donahue et 
al. (2001) and Popesso et al. (2004), the steepness 
of the X-ray to optical relations point to the fact 
that the mass-luminosity {M/ Lopt) ratio of galaxy 
clusters is not likely to be constant. This is actu- 
ally demonstrated by Popesso et al. (2005), who 
found M/L oc M"-^. Wc note that this conclu- 
sion could not be made if we had considered the 
bisector solution. Had we used the bisector, we 
would find Lx oc L^p^, which is consistent with 
a constant mass-luminosity ratio. In a future pa- 
per we plan to investigate in detail this relation 
{M/Lopt) and its possible dependence on substruc- 
ture in clusters. 

7. Summary 

We have compared a list of X-ray galaxy clus- 
ters selected from the literature (BAX) to an op- 
tical galaxy cluster catalog (NoSOCS). This com- 
parison covers the whole Northern Hemisphere for 
|6| > 30° and spans the redshift range 0.05 < z < 
0.40. The X-ray and optical centroids show ex- 
cellent agreement, with a typical X-ray— optical 
offset of < 0.50 /i-i Mpc (or < 400"). The over- 
all recovery rate of X-ray clusters in the optical 
is 81%, with the missing clusters typically poor 
and/or distant. Thus, NoSOCS efficiently recov- 
ers nearly all X-ray luminous clusters from BAX 
{Lx ^ 3.2 10^'' ergs s~^) at all redshifts below 
z ~ 0.2. In terms of richness, the recovery rate is 
90% for clusters with ^gau Z 25 out to 2: ~ 0.2, 
and for all clusters with Ngaia ^ 80 out to z 0.3. 

We employed four statistical tests to search for 
substructure using optical imaging data. These 
tests were optimized through the minimization of 
the scatter in the relations between optical and X- 
ray properties (namely Tx—^gais)- We have also 
shown the dependence of richness on the centroid 
and aperture. Wc investigated the dependence of 
substructure results on these parameters, as well 
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Tx " N^.l; (BAX coords; 57 clusters) 

Tx « (BAX coords; C i 5.0; 57 clusters) 




Fig. 18. — The connection between temperature and richness or optical luminosity. The thin lines show 
the results obtained in this work, while the thick lines exhibit results from the literature. The meaning of 
the thin lines is analogous to the previous figure, but now for the Tx^^gais relation and only considering 
clusters at 0.07 < z < 0.21. The literature results (thick lines) are given for: (i) solid line: the connection 
between Tx and Lopt,i obtained by Popesso et al. (2004); (ii) dotted line: the Tx—^gc relation obtained 
by Yee and EUingson (2003); (iii) short dashed line: the T^— Ljv/ij correlation from O'Hara et al. (2006); 
and (iv) long dashed line: the Tx~^opt,r correlation from Plionis et al. (2005). 



as on the cluster contrast, richness, redshift and 
magnitude range used for sampling galaxies. The 
substructure results shown here are based on the 
largest sample used for this purpose to date (10190 
clusters). As noted by Pinkney et al. (1996) some 
tests are much more sensitive than others. For 
instance, the /3, or symmetry, test indicates that 
35% of the clusters have substructure at the 5% 
significance level (in a radius of 1.50 Mpc and 
magnitude range of m* — 1 < nir < m* + 1), while 
the FE test finds 21%, the Lee 2D test indicates 
45% and AST only 13%. 



We have also compared richness {Ngais) and X- 
ray cluster properties (Lx and Tx)- We find that 
Tx correlates better to ISigais than Lx- We ex- 
amined the potentially dominant factors affecting 
these scaling relations. We find that the correla- 
tions are most sensitive to the aperture used for 
estimating richness. For larger radii we sample 
many more cluster galaxies, but the background 
noise also becomes much stronger. We find that 
the optimal radius is 0.50 Mpc, which is in 
good agreement to that found by Popesso et al. 
(2004). Adopting this radius for estimating Ngais, 
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we find that the cluster contrast is the dominant 
source of scatter in the Lx—^gais relation. The 
presence of substructure does not affect the scat- 
ter of this relation. On the other hand, when com- 
paring Tx to the optical parameters we find that 
substructure has a strong effect on the scatter. 

Our findings also corroborate previous results 
in the literature, although we note some issues 
that could explain some of the difficulties in per- 
forming these comparisons. These are mainly re- 
lated to the contrast of the systems, as higher con- 
trast systems generally define steeper relations. 
For instance, if we compare samples that span 
different mass regimes we may find very different 
correlations, as the contrast cuts will be very dif- 
ferent. Relations derived from richness (or Liopt) 
calculated from the X-ray centroid have smaller 
slopes than ones obtained using the optical cen- 
troid, due to the smaller galaxy number density 
seen around the X-ray centroids. The use of dif- 
ferent linear regression methods can also lead to 
inconsistent results, especially for very steep and 
noisy relations. As shown in Section 6, the bisec- 
tor solution is generally less efficient in recover- 
ing the true slopes. Literature results would likely 
show much better agreement if the same linear re- 
gression methods were used. The steeper slopes 
obtained here also suggest (as have other recent 
works) that the mass-luminosity ratio may not be 
constant. The fitting method and centroid used 
are crucial for reaching this conclusion. For the 
Lx—l^opt relation the bisector solution indicates 
Lx oc Llp^, which is close to the expectations 
from a constant mass-luminosity ratio. 

In the future we plan to use both deeper optical 
and X-ray data to investigate substructure. These 
will allow us to assess the reliability of the opti- 
cal estimates, as well as to investigate the evolu- 
tion of substructure with look-back time. We also 
plan to use apertures that scales with mass, such 
as i?5oo and i?200 to see if the structure segrega- 
tion we show here is also found with these aper- 
tures. Such virial radii are more typically used 
when measuring cluster parameters from large- 
scale simulations, and it is crucial to use the same 
techniques in comparing these with observations. 
Finally, we plan to investigate the mass-luminosity 
ratio, as well as estimate the cluster mass function 
with the NoSOCS data. 
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